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PACS 42.50.Lc - Quantum fluctuations, quantum noise, and quantum jumps 

PACS 42.50.Ct - Quantum description of interaction of light and matter; related experiments 

PACS 42 . 65 . k - Nonlinear optics 

Abstract -Non-Gaussian states represent a powerful resource for quantum information protocols 
in the continuous variables regime. Cat states, in particular, have been produced in the motional 
degree of freedom of trapped ions by controlled displacements dependent on the ionic internal 
state. An alternative method harnesses the Kerr nonlinearity naturally present in this kind of 
system. We perform detailed calculations confirming its feasibility for typical experimental condi- 
tions. Additionally, this method permits generation of all other complex non-Gaussian states with 
negative Wigner functions resulting from Kerr nonlinear interaction. Especially, superpositions of 
several coherent states are achieved at a fraction of the time necessary to produce the cat state. 
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Introduction. — Quantum information processing re- 
lying on physical observables with continuous spectra 
(continuous variables, or CV) offer new possibilities of im- 
plementation and employment when compared to their 
discrete alternatives. In particular, states with non- 
Gaussian Wigner function are essential for the speedup 
and universality of CV quantum computation [T}{3]. Ac- 
cording to the Gottesman-Knill theorem, quantum com- 
puting based only on components described by quadratic 
Hamiltonians, Gaussian inputs, and measurements on 
canonical variables can be simulated by a classical com- 
putation [J. Thus, non-Gaussian states allow to exploit 
the advantage of the quantum algorithms over the clas- 
sical ones 0. Construction of a CV universal quantum 
computer for transformations that are polynomial in those 
variables requires at least a cubic (Kerr) nonlinear opera- 
tions [J. Non-Gaussian states also find important applica- 
tions in high precision metrology [5] , novel tests of Bell-like 
inequalities [Sill], and fundamental investigations [5HTU]. 
Therefore, finding the effective and deterministic way of 
generation and manipulation of non-Gaussian states and 
interactions is crucial for further development of quantum 
based applications as well as for better understanding of 
physics. 

Position and momentum of trapped ions represent a 
good CV quantum system candidate, especially taking 
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into account the recent achievements regarding low heat- 
ing rates JT]. So far, this system is much less exploited 
than the ionic discrete variables [T^l H3] ■ In fact, mo- 
tional degrees of freedom have been mostly used as qubits, 
fulfilling the role of a quantum information 'bus' among 
the ionic internal states [M]. However, recently proof- 
of-principle quantum walks using the quadratures of the 
harmonic oscillator state were performed [T51HB] and CV 
quantum simulations were demonstrated [17] ■ 

The current ion trap technology allows an arbitrary 
quantum manipulation of the motional states |18H20j . 
In practice, however, the fidelity of operations changing 
the phonon number stepwise (e.g. blue sideband pulses) 
are comparatively low, limiting the maximum achievable 
phonon number to just a few quanta in a coherent super- 
position |21j . Alternatively, 'single-shot' methods relying 
on a single laser pulse to directly produce the desired final 
state seem to have an advantage from the experimental 
point of view. For instance, bichromatic quantum gates 
allow for high fidelities in the generation of two-qubit en- 
tangled states based on precise manipulation of the vi- 
brational degree of freedom [22H24] . as well as coherent 
superpositions involving tens of phonons |15U16j . 

In this paper we investigate a deterministic method to 
directly harness the Kerr nonlinearity present in the mo- 
tion of trapped ions and thus to produce non-Gaussian 
vibrational states |25j . The method by employing a sin- 
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gle laser pulse allows for creation of several types of 'quasi 
macroscopic' coherent state superpositions, such as the cat 
states, and all other non-Gaussian states emerging from 
Kerr nonlinear interaction. While operations capable of 
producing such states exist in principle, we believe the 
simplicity of our scheme will further push the border of 
what can be experimentally realized with high fidelities in 
CV systems. 

This paper is organized as follows. First we introduce 
the Kerr non-Gaussian state. In the next section we 
present the method of direct generation of this state on 
a vibrational mode of a trapped ion. Then the ideal Kerr 
state is compared to the result of the calculation using the 
complete Hamiltonian of the trapped ions. Wc finish the 
paper with our conclusions. 

The Kerr state. — The lowest order Hamiltonian 
in powers of the harmonic oscillator annihilation a and 
creation a) operators capable of producing non-Gaussian 
states describes a Kerr nonlinear medium. The evolution 
of an initial coherent state \a) subject to this Hamiltonian 
gives rise to the so-called 'Kerr states', 
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Fig. 1: Evolution of the ideal Kerr state, (a) Averages and 
variances of the quantum oscillator position X and momentum 
Y operators as a function of the effective time r [26] . The insets 
below show a few Wigner functions for specific times [27]: (b) 
Initial spreading in phase space, r = n/5; (c) Superposition of 
three coherent states, r = 2ir/3; (d) Cat state, r = n. 



They are parametrized in terms of only one quantity, the 
effective time r [55] . In general, the Kerr states are highly 
nonclassical and for certain values of t their Wigner func- 
tion assumes negative values [37]. Their most prominent 
examples are the cat states, formed by superpositions of 
two coherent states shifted in phase by tt. For instance, 
the cat state e~ ln ^\ia) + e I7r / 4 | — ia) results from Eq. (fTJ 
taken for t = tt. After this time the evolution effectively 
reverses to reach the original coherent state at r = 2tt. 
Cat states were used in studies of decoherence mechanisms 
and the quantum-classical boundary [8] [9]. The 'large 
cats' for which the two components \ia) and | — ia) are 
nearly orthogonal (a> 1.5) also find applications in quan- 
tum information processing [55] US]- Efforts were made 
to produce such states using the Kerr nonlinearity in op- 
tical fibers [30]. However, optical nonlinearities are far 
too small to enter the highly nonlinear regime and thus 
to obtain nonclassical Wigner functions before dissipation 
effects destroy the coherences. Superpositions of multiple 
coherent states are obtained for r = 7r/2n, where n is the 
number of coherent states in the superposition [3TU32] . For 
instance, at r e s = n/2 one achieves the superposition of 
four of them (the compass state) and for r e g = ir/3 six co- 
herent states are superposed |33j . The maximum number 
of participating coherent states relies solely on the physical 
possibility of distinguishing them, which depends on the 
amplitude of the initial coherent state. Larger values of a 
will allow more complex superpositions to be constructed, 
but will have support on a larger number of Fock states as 
well. Such examples of Kerr states were not experimen- 
tally investigated so far. 



Time evolution of mean values and variances of the posi- 
tion X = a + a^ and momentum Y = — i(a — a?) operators 
in Kerr medium is depicted in fig. [T^, ■ The evolu- 

tion spreads the initial coherent state (A 2 X = A 2 Y = 1) 
around the origin of phase space, soon producing negative 
values in the Wigner function. Negative values with am- 
plitude around 10% of the maximum positive value are al- 
ready observed for t ps 7t/20. In fig. [Tp they reach nearly 
the same magnitude in sub-Planck areas of phase space 
for t ps 7r/5. Since the Kerr interaction does not cou- 
ple different Fock states, but instead only dephases them 
independently, the Wigner function cannot occupy higher 
phonon numbers than originally present. For intermediate 
times, fig. [Tp shows one possible superposition of multiple 
coherent states, forming in this case a triangular pattern. 
The cat state is obtained at the turning point of the evolu- 
tion (fig.[JJi), after which the states produced appear once 
more in reversed order. 

Generation of motional non-Gaussian Kerr 
states. — Let us consider a single ion trapped in the ef- 
fective harmonic potential of a Paul trap interacting with 
a resonant laser field. By proper choice of geometry, only 
one of the three existing vibrational modes can be made 
relevant to the dynamics. In this case the interaction 
Hamiltonian reads 1351 



H = h 



n 



(2) 



where f2 is the Rabi frequency, er + is the electronic state 
rising operator, u is the quantum oscillator harmonic fre- 
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quency, and r\ is the Lamb-Dicke parameter. We will fur- 
ther assume the Lamb-Dicke regime and thus (a^a)rj < 1. 
It is justified in this case to expand the exponentials of 
eq. ([2]) in rj and disregard the higher order terms. The 
Kerr effect is described by the term (a'a) 2 . It dcphascs 
each Fock state proportionally to its eigenvalue (self-phase 
modulation) [25] . Keeping only the terms resonant to the 
atomic transition (carrier) up to the forth order in 77, we 
obtain in the interaction picture 



H int = ft 
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The first three terms describe the corrected internal tran- 
sition Rabi frequency, and do not influence the motional 
state. The terms proportional to a^a represent a rota- 
tion in phase space and can be omitted by considering an- 
other rotating frame. Therefore, the lowest order term to 
nontrivially influence the motional dynamics is indeed the 
Kerr self-phase modulation. Higher order terms, although 
contributing as well, are much smaller by the factor of or- 
der of rf . Similarly, off- resonant terms will also hinder 
the ideal dynamics. All these effects are disregarded in 
the present qualitative discussion, but are fully taken into 
account in the numerical calculation presented in the next 
section. 

The inclusion of additional vibrational modes brings 
other nonlinear effects such as cross-phase modula- 
tion |36j . Here we will disregard those interesting but more 
complicated effects by supposing that the Lamb-Dicke pa- 
rameters of such modes are much smaller than for the main 
mode considered. If more ions participate in the dynam- 
ics, all the vibrational modes in a given direction must be 
included for consistency, in which case cross-phase modu- 
lation is unavoidable. 

The Hamiltonian of eq. ([3]) couples the internal and 
the external degrees of freedom and results in an inter- 
nal state-dependent evolution of the motion. To remove 
this undesired effect, the atomic internal state must be 
initially prepared in an eigenstate of a x . Usual 7r/2 Rabi 
pulse with 7r/2 phase with respect to the main Kerr pulse 
accomplishes the atomic state superposition on the time 
scale of f2 _1 , i.e. much faster than the Kerr evolution 
time scale r e g of eq. (|5|). Thus, it has a negligible effect on 
the ion motional state. With this procedure, the internal 
atomic state is separable from the motional state. The 
evolution operator acting on the vibrational mode then 
equals 

U(t) = cxp jiifijjVVt j . (4) 

From this expression we identify the effective value of the 
Kerr evolution parameter 
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can be prepared in a number of ways [T2JQ3] ■ Once pre- 
pared, any state can be reconstructed by a standard tech- 
nique using displacements and parity measurements |37j . 

A special feature of the method under discussion is the 
independence between the cat state 'size' (with respect to 
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and the time necessary to produce it. Since 



The highly non-Gaussian Kerr states require a coherent 
state as the starting point of the nonlinear evolution. It 



the Kerr evolution recurrence time does not depend on the 
initial state, a 'small cat' and a 'large cat' are generated 
by the same light pulse. This is in sharp contrast to ap- 
proximative methods where the desired final superposition 
is composed of several laser pulses each rising the phonon 
number by one unity |18| . In the scheme we investigate 
here, the maximum cat size is only limited by the trap 
anharmonicities and decoherence (as all other methods ul- 
timately are). Furthermore, no additional truncation at a 
maximum phonon occupation number is necessary. 

A simple estimate using eq. (|5|) shows the feasibility 
of the scheme within the typical experimental conditions. 
We take ui = 2ir x 3 MHz as the quantum oscillator fre- 
quency. The Rabi frequency 51 — 2tt x 200 kHz would 
result in off-resonant excitation on the order of 0.1%. For 
these parameters, ?y = 0.1 would allow achieving the cat 
state in i cat ~ 100 ms. This value lies on the boundary of 
the expected motional coherence time, typically less than 
1 ms for small traps and above 100 ms for larger traps. 
However, because of the strong dependence of t ca t on the 
Lamb-Dicke parameter, 77 = 0.2 would already result in 
a pulse duration of only i cat ~ 6 ms. In any case, other 
complex Kerr states would be accessible for even shorter 
pulses, as previously discussed, lying well within the ex- 
pected motional coherence time. The Lamb-Dicke param- 
eter can be manipulated to a certain extent by changing 
the oscillator frequency as well as by small adjustments in 
the angle between the oscillator axis and the laser propa- 
gation vector (or two lasers, in the case of Raman excita- 
tion) I2S]. 

Fidelity of the Kerr state produced by the com- 
plete Hamiltonian. — The complete internal and ex- 
ternal dynamics of the trapped ion using the full Hamil- 
tonian of eq. ([2]) was computed. The effects of the higher- 
order and off-resonant terms were included in the calcula- 
tion. As the initial condition for the ion motion a coher- 
ent state of amplitude a = 2 was taken, corresponding to 
mean phonon number (a^a) = 4. The oscillator frequency 
ui = 2tt x 3 MHz was assumed. 

For the numerical calculation the Fock basis was trun- 
cated at \n = 21) in order to obtain the probability smaller 
than 10" 19 for the population of the higher Fock states. 
The initial electronic state (\g) + |e))/ y/2 was chosen as the 
eigenstate of a x (where \g) and |e) stand for the ground 
and the excited internal state, respectively). The evolu- 
tion of the vibrational state was calculated by applying the 
evolution operator based on the Hamiltonian of eq. (J2J in 
steps of dt = 10" 5 ms or smaller. The states obtained 
in this manner were compared to the ideal Kerr state of 
eq. ((1} using the fidelity defined as the square of their 
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Fig. 2: (a) Evolution of the first moments of the oscillator 
position X and momentum Y, starting from the coherent state 
\a = 2), for r) = 0.2, fi = 2tt x 100 kHz, and oj = 2n x 
3 MHz. Compare to fig. [1] (line colors and styles follow the 
same convention). Physical time is displayed on the lower axis, 
whilst the effective Kerr time appears on the upper axis. The 
fidelity of the motional state evolution as compared to the ideal 
Kerr state of eq. Q]is depicted by the gray curve (right axis), 
(b) Internal state population, (c) Phonon population. 



scalar product. 

An example of pulse operation as a function of time for 
i] = 0.2 is shown in fig. [21 For this particular choice of pa- 
rameters, the cat state was obtained in 13 ms interaction 
time with 98.5% fidelity. More complex states were gen- 
erated faster with higher fidelity. The fast ripples on the 
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Rabi frequency (fl/co) 

Fig. 3: Cat state fidelity as a function of Q, for three different 
values of rj (full symbols). The times needed to produce the 
cat state are displayed on the right axis (open symbols). The 
dashed lines result from the approximated effective evolution 
time of eq. [S] Each point results from a numerical calculation 
of the quantum evolution. 



internal state population (fig. ^jp) reach maximum ampli- 
tudes of 1% (although typically 0.1%), whilst the phonon 
state populations (fig. [2}:) fluctuates within 0.1%. Those 
off-resonant excitations, however, do not accumulate dra- 
matically over time. Those results confirm that the dy- 
namics using higher-order terms of the Hamiltonian docs 
not hinder in itself the attainability of high fidelities. 

The cat state fidelity as a function of Rabi frequency for 
three values of r] is depicted on fig. |3J The smaller values 
of rj, the higher fidelities for the cat state. The physical 
time necessary for the formation of the cat state is shown 
on the right axis. In principle, fidelity of the state can be 
arbitrarily close to unity although the time required for 
its creation increases. In practice, decoherence will limit 
the maximum operation time. Thus, the optimum pulse 
duration will depend on the particular experimental con- 
ditions. The tradeoff between increasing Rabi frequency 
and decreasing the ideal fidelity results in shorter pulse 
duration. In fact, the operation could be faster even for 
low values of rj if the Rabi frequency would be increased 
to values comparable to the oscillator frequency. It is im- 
portant to note however, that eq. ([5]) is no longer valid in 
the limit of small rj and large due to saturation effects 
and off-resonant excitations. 

The more realistic experimental situation where the cat 
state generation has to be completed given a fixed inter- 
action time was also examined. The cat state fidelity as a 
function of the Lamb-Dicke parameter rj for a 10 ms pulse 
duration is shown in fig. [4] Rabi frequencies are indicated 
on the right axis of the figure. The optimum cat state 
fidelity of 99% was found for rj « 0.15 for our choice of 
parameters. The smaller r\ is, the lower the fidelity be- 
comes due to the large values of fi which drive undesired 
transitions. As expected, the higher order terms in eq. ^ 
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Fig. 4: Cat state fidelity as a function of the Lamb-Dicke pa- 
rameter r\ for a fixed cat state creation time of 10 ms (black 
circles) . The required Rabi frequencies relative to the oscillator 
frequency (Q./ui) are indicated by the red squares. Each point 
results from a numerical calculation of the quantum evolution. 

contribute significantly for larger r\ and hinder the ideal 
Kerr dynamics. 

Although the Kerr Hamiltonian allows one to generate 
cither 'small' or 'big' cats using the same interaction time, 
bigger cats will decohere faster, since coupling to a thermal 
reservoir results in loss of visibility of the interference pat- 
tern in phase space (with exponential dependence on the 
cat size) [33] ■ However, the cat state is the final product of 
the Kerr interaction, and as such takes the longest time to 
be created. Other interesting macroscopic quantum super- 
positions could therefore coherently occupy larger areas in 
phase space albeit for shorter durations. 

Comparison between the quantum evolution of some of 
the states presented in fig. 0] and the ideal Kerr states is 
presented in fig. [3J For the optimum i], a superposition of 
four coherent states is produced with 99.6% fidelity, whilst 
three superposed coherent states can be obtained with 
99.3% fidelity. The first negativities of the Wigner func- 
tion are observed already after 500 fj,s interaction time, 
and dramatic effects in sub-Planck regions of phase space 
appear after 2 ms with more than 99.9% fidelity. We stress 
that these nearly ideal fidelities indicate that the actual 
experimental limitations will not be a consequence of the 
form of the Hamiltonian, but rather decoherence and im- 
perfections which will inherently depend on specific details 
of the experimental setup. For instance, larger cats can 
be expected to experience stronger decoherence [35] , effec- 
tively requiring faster nonlinear dynamics. 

From the perspective of the interaction Hamiltonian, 
the upper limit on the fidelity for a fixed cat state cre- 
ation time results either from off-resonant excitation in- 
volving the phonon sidebands or from higher-order terms 
in the Hamiltonian. Off-resonant excitations cause noisier 
evolution for smaller ?y, while for larger values of rj the 
evolution becomes smoother (fig. [5]). These excitations 



Fig. 5: Fidelity of the motional state evolution as compared to 
the ideal Kerr state. A fixed time of 10 ms (corresponding to 
t — 7r) is imposed for the formation of the cat state (fig. 2J. 

can be thus avoided by increasing the oscillator frequency. 
Taking all necessary compensations into account, our cal- 
culations indicate that a higher oscillator frequency has a 
positive effect on the state fidelity. 

Conclusion. — We have presented a detailed inves- 
tigation of a practical method of involving Kerr non- 
Gaussian state generation in the vibrational mode of 
trapped ions [55] . The method employs only one laser 
pulse, directly connecting an initial coherent state to the 
desired highly nonclassical Kerr state. Since only a reso- 
nant carrier pulse is needed to harness the natural Kerr 
nonlinearity of trapped ions, this scheme can be performed 
even in very simple setups, such as the stylus trap [40] . 
An optimized setup would be desirable, however, where 
the trap size should be large enough to avoid motional 
heating and decoherence at the same time that a high 
oscillator frequency would hinder off-resonant excitations 
even for fast dynamics. 

Special attention has been paid to the cat states, for 
their importance in quantum information applications. 
However, other complex quantum states such as super- 
positions of several coherent states can also be generated 
in the process. Highly non-Gaussian states with negative 
values on the Wigner function are obtained for very short 
interaction times. 

The scheme produces large superpositions in phase 
space using the same time resources necessary for small 
superpositions, without the need for artificially truncat- 
ing the population in the Fock state basis. 
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